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Nonperturbative Results for Yang-Mills Theories 
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Some non perturbative aspects of the pure SU(3) Yang-Mills theory are investigated assuming a 
specific form of the beta function, based on a recent modification by Ryttov and Sannino of the known 
one for supersymmetric gauge theories. The characteristic feature is a pole at a particular value of 
the coupling constant, g. First it is noted, using dimensional analysis, that physical quantities behave 
smoothly as one travels from one side of the pole to the other. Then it is argued that the form of 
the integrated beta function g(j-i), where p is the mass scale, determines the mass gap of the theory. 
Assuming the usual QCD value one finds it to be 1.67 GeV, which is in surprisingly good agreement 
with a quenched lattice calculation. A similar calculation is made for the supersymmetric Yang-Mills 
theory where the corresponding beta function is considered to be exact. 
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I. INTRODUCTION 

The exact beta function of a gauge theory is gener- 
ally considered to contain many non-perturbative secrets 
of the gauge theory behavior. Unfortunately it seems 
to be only computable analytically in perturbation the- 
ory. Physically it is related to the trace anomaly, or the 
non-zero value of the divergence of the scale (dilation) 
current. We are specializing to massless theories here so 
it represents a violation of the classical result. 

Another object in gauge theories with massless 
fermions, the divergence of the U(l) axial vector current, 
should be zero in the classical limit but is known not to 
vanish at the quantum level (axial anomaly). Most inter- 
estingly it has the property that there are no corrections 
[1 1 beyond the one loop level so it can be considered to 
be known exactly. 

Now, in supersymmetric gauge theories, it was found 
[2 1 that both the trace anomaly and the axial anomaly 
(actually together with the divergence of a special su- 
perconformal current) belong to the same chiral multi- 
plet and hence should be somehow related to each other. 
There was a lot of discussion about the meaning of this 
feature - for example is there a contradiction between an 
exact one loop result for the axial anomaly and a result 
containing all orders of coupling constant for the trace 
anomaly - and finally it was realized that there could 
be compatibility, with the exact beta function being de- 
termined from the axial anomaly. This was shown both 
for the supersymmetric Yang-Mills theory |3 1 and for the 
supersymmetric theory containing also fermions not be- 
longing to the adjoint representation |4|. A characteristic 
feature of these exact beta functions is a pole at a par- 
ticular value of coupling constant, g in the beta function 
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Recently, Ryttov and Sannino [5] conjectured such all 
order beta functions for ordinary SU(N) (non supersym- 
metric) gauge theories both with and without fermions 
belonging to arbitrary representations based on analogy 
to the supersymmetric case. The generalization, by one 
of the authors, for the SO, and Sp gauge groups ap- 
peared in ] 6 1 and for chiral gauge theories in [7\. These 
beta functions were found to satisfy known consistency 
conditions at second order and to work well in many 
interesting applications to working technicolor models 
reviewed in[7]. These beta functions also feature a pole 
at a particular value of g. 

In the present note, we study some non-pertubative 
consequences of the conjectured beta function for the 
simple SU(N) Yang-Mills theory; already |5J the conjec- 
tured beta function had been seen to give a reasonable 
picture in the asymptotically free perturbative region. In 
section II we use dimensional analysis to investigate the 
running with scale of physical quantities with various 
engineering dimensions. This involves g(j.i) but it was 
noted that physical quantities change smoothly as g goes 
through the pole value. 

In section III, we investigate the integration of the 
defining equation for the beta function which yields an 
explicit expression for g(pi)- We consider all Yang-Mills 
theories in the sense that a full range of values for the 
coupling constant at a reference mass are considered. 
An amusing feature is seen to arise: the solutions for 
g(jj) do not allow /.( to be lower than a certain value. 
We interpret this as the measure of the mass gap for the 
Yang-Mills theory. For the usual value of the QCD cou- 
pling constant, with N=3, our predicted value of 1.67 
GeV is seen to be in good agreement with a valence lat- 
tice calculation. We also compute the mass gap for the 
supersymmetric SU(3) and SU(N) Yang-Mills theories. 
In the latter case, the starting beta function does not re- 
quire any conjecture. 

Section IV contains a brief summary and discussion. 
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II. GENERAL SCALING RESULTS FOR YANG - MILLS 
THEORIES 

Consider a renormalization group invariant quantity 
H of mass dimension D in any theory with a single field 
and a single coupling constant g. Define p. to be the 
renormalization scale and gip) to be the value assumed 
by the coupling constant at that scale. Since H must 
be a measurable physical quantity dimensional analysis 
implies that it has the form, 



H = p D <H[gip)] 



(1) 



where "H [gip )] is some function of the coupling con- 
stant. Of course, for any given theory there are many 
different interesting quantities H. One may, as usual, in- 
troduce a characteristic invariant scale, A for the theory by 
defining a particular H to be A D . By differentiating both 
sides of Eq.tfTl with respect to log(f(), the natural log, and 
recognizing that the left hand side is independent of p, 
we obtain the main equation: 



, with p{g) 



dlogiy.) 



(2) 



By integrating this equation one immediately finds: 

J 2a 



Hence if the exact beta function were known as a func- 
tion of g, one could express any renormalization group 
invariant quantity H as a known function of g up to an 
arbitrary overall constant. 

Now, as discussed in the preceding section, a conjec- 
tured all orders beta function for the ordinary Yang-Mills 
theory based on the known supersymmetric Yang-Mills 
theory all orders beta function was recently introduced in 
[5 1 by Ryttov and Sannino. It was found to be consistent 
with other non perturbative approaches to the ordinary 
(non supersymmetric) gauge theories and hence to be a 
reasonable model for further investigation. 

This RS all orders beta function ansatz for the SU(N) 
Yang-Mills theory reads: 



there is a different numerical constant for each choice 
of H. These numerical constants might be approximated 
by using perturbation theory at a large value of p, for 
example. 

For definiteness we list some possible interesting 
choices for H. 

i. The gluon condensate or vacuum expectation 
value of the trace of the energy momentum ten- 
sor. This has the engineering dimension D = 4. 

ii. The glueball squared masses (D = 2). We presume 
that there is a spectrum of glueballs with different 
spin-parities and masses. 

iii. The glueball-glueball scattering cross sections (D 
= -2). 

iv. The coefficients, a„ in the expansion of the partial 
wave amplitudes for glueball scattering, Ha n s n in a 
region of analyticity (Here D= -2n). 

To get an idea of the dependences of various phys- 
ical quantities on the coupling constant g we plot the 
characteristic factor, 



>YM 



with a 
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Using Eq.Q in Eq.Q then yields: 



and b 



17 N 
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where the overall constant is defined by 

const = fiogo ^ exp ' 



2a; 



(5) 



(6) 



Here the subscript zero denotes the value corresponding 
to the lower limit of integration in Eq. ||3}. Of course, 
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exp 
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(7) 



in Eq.(|5j for the case D = 4 in Fig. [I] (a) one and for the 
case D = -2 in Fig. [T] (b) . For definiteness the choice N 
equal to 3 is made. It is seen that each of these curves 
displays an extremum at the same value of g. 
Differentiating shows that this value is, 



= 1/b = 



88n 2 
51 



(4.127) 2 



(8) 



Furthermore, the second derivative of Fig) at this point 
is simply — 2DF(g) lag 4 ; this means that dynamical quan- 
tities with positive dimension D will have a maximum at 
this point while those with negative D will have a mini- 
mum at the same point. We also note the interesting fact 
that for a given value of the physical quantity F, there are 
two different values of g. The situation will be further 
explored by looking at the beta function. 



III. BETA FUNCTION AND MASS GAP 

A graph of the N =3 Yang-Mills conjectured beta func- 
tion in Eq.Q is shown in Fig. 2(a). From its shape one 
concludes mat the origin, g =0 is an ultraviolet stable 
fixed point and that there is clearly a pole at g = b~ 1 ^ 2 . 
However, we observe from the Figs. [11(a) and (b) that the 
physical quantities, H remain smooth as one goes from 
one side of the pole to the other. 

In order to discuss the coupling constant in a gauge 
theory, one chooses a reference value jUo of the energy 
scale and specifies (from experiment) the value of the 
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(a) Case D = 4 




(b) Case D = -2 
FIG. 1: F(g) as function of g for N = 3. 



(a) [>YM(g) a s function of g for N = 3. 




(b) log(fi/fio) as function of the coupling g 

FIG. 2: Plots of the YM beta function, panel (a), and of the solu- 
tion for the running of the scale (b) as function of the coupling 
constant g with g = 1. 



running coupling constant, g(fio) at this point. In the 
case of QCD, it is convenient to choose [Iq - niz ~ 91.19 
GeV. For illustration we will also choose this value of po 
as our reference scale. The value of the QCD coupling 
constant at this scale is measured to be about 1.228. In 
the present pure Yang-Mills case we are dealing with a 
hypothetical model so we are free to choose any value 
for g(po). In fact it is better to allow g(po) to range over 
all possible values, which corresponds to describing all 
pure Yang-Mills N=3 theories. 

Clearly, the first step is finding out how the coupling 
constant runs with the scale p. Integrating the second of 
Eq.S yields the following relation between g = g(p) and 
log(p) corresponding to any choice of go = gi^o)'- 



" ,8l iK ?-= 



« \go 



(9) 



It is convenient to plot in Fig. 2(b) log (fi/fio) as a func- 
tion of g with the choice go = 1 . In the part of this plot 
to the left of g([i) = b~ 05 « 4.127, g(p) decreases as log(/.() 
increases. That is the expected asymptotically free be- 
havior. On the other hand, to the right of this point, 
g(p) smoothly starts rising with increasing log(/.(). There 
is no discontinuity at the pole. Clearly, this behavior is 
the same as that shown for physical quantities in panels 



(a) and (b) of Figs. [T] The existence of a smoothly con- 
nected different phase for the theory is intriguing. For 
the present, however, we will concentrate on the asymp- 
totically free region. 

The different N = 3 pure Yang-Mills theory which is 
defined by go = 3 yields the running coupling constant 
plotted in Fig. 3(a). The overall picture is substantially 
the same. In particular, the asymptotically free region 
still corresponds to g(p) < b~ 05 . Similarly the theory 
characterized by go =0.5 is seen from the plot in Fig. 3(b) 
to still have the same range in g(p) for the asymptotically 
free region, but with rather different p values. 

Figures 2(b), 3(a) and 3(b) each illustrate that there is 
a particular (but different in each case) value of log(/i) 
below which the curve does not extend. This is peculiar 
since it would imply that the running coupling constant 
could not be measured experimentally in that region. The 
only way in which this might be consistent is if the disal- 
lowed region would be lower than M, twice the mass of 
the lightest glueball state in the theory. Then it would be 
below threshold and not accessible to experiment. That 
seems like a plausible determination of the mass gap of 
the theory. It is clear that the value of such a mass gap is 
determined by the ordinate of the point where log(f() is 
minimum, i.e. where g = b~ 05 . Specifically, 



3 



CP 3 - Origins 



Particle Physics & Origin of Mass 




(a) £o = 3 




(b) go = 0.5 

FIG. 3: Solution for running scale vs. coupling constant g. 



1 f 1 

log(M) = log^o) + Y a \ h ~~ 



OJ 



-log(g b°% (10) 



where go is the assumed value of the coupling constant 
at fig = 91.19 GeV. 

Note that the so defined mass gap depends on the par- 
ticular N=3 Yang-Mills theory we are considering via its 
dependence on go- In Figure [4] (a) twice the mass gap, 
M in GeV is plotted for g ranging from 1 to 10. For 
orientation, the mass gap for go slightly larger than one 
is seen to be in the GeV range. The mass gap reaches a 
maximum when go takes the pole value, 1/b 05 and de- 
clines somewhat for larger values of go- It seems reason- 
able that within the asymptotically free regime, the mass 
gap increases with the strength of the coupling constant, 
go. For go < 1, the mass gap starts to decline extremely 
rapidly, as may be seen from Fig. |4](b). 

To further test the reasonableness of our interpretation, 
we may try to predict the mass of the lightest glueball in 
this model. We would like the result to be similar to the 
usually expected value (about 1.5 GeV) when we adopt 
the QCD value for the coupling constant, go - 1.228, 
mentioned above. This would correspond to neglecting 
the effects of quark fields on the glueball mass. Since it 
is seen that there is a rapid dependence of M on go we 
give in Fig. [5] a blow up of the prediction for M/2 in the 
region of go near the experimental one. This yields for 



the lightest glueball mass, 

M/2 « 1.67 GeV, 



(11) 



which does seem reasonable for an a priori prediction. 
We may compare this value with the result of a lattice 
QCD calculation employing a valence assumption |8J. 
That calculation gave the result for the + glueball which 
is found to be lightest: 



m(0 + ) = 1648 ± 58 MeV. 



(12) 



This embarassingly accurate agreement gives us at least 
some confidence in the correctness of the interpretation 
of the mass gap and the validity of the RS conjectured 
beta function. Other lattice results are also in agreement 
with the value quoted above ||9llT0l . To sum up, the pole 
in the beta function does not produce any singularity in 
the theory but seems to be the feature which generates 
the mass gap. It would be interesting to investigate the 
dependence of this result on the choice of renormaliza- 
tion scheme. 

If the conjectured pure Yang-Mills beta function could 
be shown to be the same as the exact one, it would 
amount to showing that the theory has the kind of mass 
gap defined above. 

Actually, it is easy to see that this mass gap mechanism 
does not, in general, require the RS conjecture if one 




(a) M (twice the mass gap) dependence on 1 < gg < 10 




(b) M (twice the mass gap) dependence on go < 1 

FIG. 4: M (twice the mass gap) dependence on gg for j.ig 
91.19 GeV and g = b~ os . 
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FIG. 5: Mass gap dependence on 1.1 < g < 1.3 for jj. 
91.19 GeV and g = 2r°- 5 . 




FIG. 6: Solution of the running of the scale as function of the 
coupling constant g for the Super Yang-Mills theory. 



goes beyond the ordinary Yang-Mills model. That is 
because the supersymmetric Yang-Mills theory, which 
is of course the starting point of the RS conjecture, is 
known to possess an exact beta function of the same 
form as Eq.Q but with somewhat different values of a 
and b, namely: 



the ordinate at the minimum point) is shown in Fig. [7]( 
It shows that the mass gap has a maximum as function 
of the number of colors. 

When considering the large N behavior of gauge the- 
ories it is often desired to make an extrapolation where 
large N is taken in such a way that the quantity, 



a = 



3N_ 

(4ti) 2 



N 

8n 2 ' 



(13) 



The numerical results are qualitatively similar if one 
assumes a similar value of the coupling constant. For 
example if we choose to specify the super Yang-Mills 
theory by taking a similar coupling constant at the scale, 
jj-o = niz, we would get for the choice go = 1.0, the curve 
shown in Fig. [6} which describes the running of g(f.i) with 
respect to log(p) and can be compared with the curve in 
Fig. 2(b). 

We can compare the mass gap, M'/2 in the supersym- 
metric case by using Eq.([T0| and making the choice go = 
1.228 (and also N=3). This would give for the mass of 
the lightest supersymmetric multiplet, 



(15) 



is held fixed. With such an extrapolation there will 
be no dependence (in either the supersymmetric Yang- 
Mills or the pure Yang-Mills with RS type beta function 
cases) on N. This may be immediately seen by writing, 
from Eq.ljijl, a = aN and b = bN so that a and b are 
independent of N. Then, as an example, Eq.jlO) for the 
mass gap becomes 



log 



M 



b (1 
S\2 



l0 ^ N, ' 2) )-2 S ^N' ( ' 6) 



which only involves N via the fixed combination g^N. 



M'/2*0.49 GeV. 



(14) 



Clearly, the supersymmetric and non supersymmetric 
Yang-Mills model results seem to be only qualitatively 
similar, the lightest supersymmetric particle having a 
mass about 1/3 that of the non-supersymmetric model 
glueball (assuming the same coupling constant). 

Of course, another interesting aspect to explore for 
gauge theories is their behavior as the number of colors, 
N gets large. Taking the supersymmetric gauge theory 
as an example and fixing, for the sake of definiteness, 
the gauge coupling constant as 1.228 at j.io = 91.19 GeV, 



we find the running of g(f.i) by substituting Eq.l 13 1 into 
Eq.([9j. Figure [71 (a) show log(/j) plotted against g(j.i) 
for respectively N=3 (blue solid line), N=15 (red dashed 
line) and N=100 (black dot-dashed line). It is seen that 
the asymptotically free (left) region in g shrinks as N 
increases. Furthermore, the mass gap (corresponding to 



IV. SUMMARY AND DISCUSSION 

First, an exploration of the scale dependence of phys- 
ical quantities (with engineering dimension, D) was 
made using dimensional analysis and a specific form, 
conjectured by Ryttov and Sannino, of the beta function 
for the pure Yang-Mills theory. It was noted that even 
though the beta function had a pole (inherited from the 
known form for supersymmetric gauge theories which 
stimulated the conjecture) physical quantities remained 
smooth as the pole value of the coupling constant, g was 
crossed. 

For a more detailed understanding of how the cou- 
pling constant runs, the integration of the beta function 
was next carried out for the complete set of SU(3) Yang- 
Mills theories, i.e. those corresponding to any choice 
of coupling constant, go at a convenient reference scale, 
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Yang-Mills. 



FIG. 7: In the panel (a) we plot the solution of the running of the scale as function of the coupling constant g for the Super 
Yang-Mills theory for N = 3, solid blue line; N = 15, red-dashed line and N = 100, black-dot-dashed line. We have kept fixed 
go = 1-228 at the scale fi = 91.19 GeV. In the panel (b) we plot the mass gap M'/2 as function of the number of colors N. 



f.io = niz- Then g(y) was also seen to have a smooth 
behavior at the pole value of g. Most interesting is that 
(due to the existence of the pole) the curve of g(f.i) pre- 
dicts a numerical value for the mass gap of the Yang- 
Mills theory, i.e. the mass of the lightest glueball. The 
predicted value, 1.67 GeV, seems rather close to the one 
obtained from a lattice treatment of QCD in the valence 
or quenched approximation. 

It is noted that if the reference value of the coupling 
constant, go decreases below about 1 (the experimental 
value is 1.228) the mass gap drops very quickly. 

A similar treatment was carried out for the supersym- 
metric SU(3) Yang-Mills theory and the mass gap was 
calculated to be about 0.49 GeV. In this case, the form of 
the beta function is quite similar but is based on a known 
rather than a conjectured beta function. 



We plan to next investigate the extent to which this 
work can be carried out for non-supersymmetric gauge 
theories containing fermions. Evidently there are many 
interesting questions which remain. 
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